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Heavy quark production has been very well studied over the last years both theoretically and 
experimentally. Theory has been used to study heavy quark production in ep collisions at HERA, 
in pp collisions at Tevatron and at RHIC, in pA and dA collisions at RHIC and in AA collisions 
at CERN-SPS and at RHIC. However, to the best of our knowledge, heavy quark production in 
eA has received almost no attention. With the possible construction of a high energy electron-ion 
collider, updated estimates of heavy quark production are needed. We address the subject from 
the perspective of saturation physics and compute the heavy quark production cross section with 
the dipole model. We isolate shadowing and non-linear effects, showing their impact on the charm 
structure function and on the transverse momentum spectrum. 
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Keywords: Quantum Chromodynamics, Heavy Quark Production, Saturation effects. 



I. INTRODUCTION 

The construction of a high energy Electron Ion Collider (EIC) was proposed in 2005 fT] (See also 0]). During 
the subsequent years, several predictions for the inclusive and difFractive observables were made, especially in the 
context of saturation physics. One of the advantages of this new machine is that it will be possible to reach values 
of the saturation scale, Qs, which are larger than those reached at HERA. A large saturation scale is crucial for the 
observation of most of the saturation effects. In particular, the collider environment is ideal for studying semi-inclusive 
and exclusive processes. In previous works, [3t3|j '^^ made predictions for the inclusive nuclear structure function F2 
and Fl as well for the diffractive observables. In these works we have concluded that the nuclear structure function 
is reduced up to 50% with respect to case where saturation effects are not taken into account. We made estimates 
for the ratio between the nuclear difFractive and total cross sections and predicted that about 30% of the events at an 
EIC will be difFractive. We have also investigated the dependence on the /3 and xp variables of the nuclear diffractive 
structure function xpF^^K We showed that xpF^f^ becomes very flat in j3 and xp when we increase the atomic 

number, A, and we found the same behavior for the ratio R = P2A1 1 ^'2A2 ^^"^ different nuclei. Concerning the 
exclusive vector meson production off nuclei, we showed that the coherent process (when the nucleus remains intact 
after the collision) for vector meson production will be much more important than the incoherent one (when the 
nucleus breaks up after the collision). 

In this paper we continue our study of quantities that could be measured in an electron ion collider and calculate 
the cross section of heavy quark production using the dipole approach and a nuclear saturation model based on the 
physics of the Color Glass Condensate (CGC) (For related studies see Refs. 8, 9]). The main input of our calculation 
is the dipole-nucleus cross section, fJdA{x,r), which is determined by the QCD dynamics at small x. In the eikonal 
approximation it is given by: 

aaA{x,r)=2 j <fbM^{x,r,b) (1) 

where M^{x^ f:t>) is the forward dipole-target scattering amplitude for a dipole with size r and impact parameter b 
which encodes all the information about the hadronic scattering, and thus about the non-linear and quantum effects 
in the hadron wave function (see e.g. [l^). It can be obtained by solving the BK (JIMWLK) evolution equation 
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FIG. 1: (Color online) Nuclear charm structure functions for A = Pb and different values of considering different models 
for the dipole-nucleus cross section. 



in the rapidity Y = ln(l/2:) ^JTs^. Many groups have studied the numerical solutions of the BK equation, but 
several improvements are still necessary before using the solution in the calculation of observables. In particular, one 
needs to include the next-to-leading order corrections into the evolution equation and perform a global analysis of 
all small x data. It is a program in progress (for recent results see fisj). In the meantime it is necessary to use 
phenomenological models for Af which capture the most essential properties of the solution. Following ^] we will use 
in our calculations the model proposed in Ref. [iBl, which describes the current experimental data on the nuclear 
structure function as well as includes the impact parameter dependence in the dipole nucleus cross section. In this 
model the forward dipole-nucleus amplitude is given by 



Af^{x,r,b) = 1-exp 



--adpix,r^)TA{b) 



(2) 



where adp is the dipole-proton cross section and TA{b) is the nuclear profile function, which is obtained from a 3- 
parameter Fermi distribution for the nuclear density normalized to A (for details see, e.g., Ref. Q). The above 
equation, based on the Glauber-Gribov formalism [13], sums up all the multiple elastic rescattering diagrams of the 
qq pair and is justified for large coherence length, where the transverse separation r of partons in the multiparton 
Fock state of the photon becomes a conserved quantity, i. e. the size of the pair r becomes eigenvalue of the scattering 
matrix. It is important to emphasize that for very small values of x, other diagrams beyond the multiple Pomeron 
exchange considered here should contribute {e.g. Pomeron loops) and a more general approach for the high density 
(saturation) regime must be considered. However, we believe that the present approach allows us to estimate the 
magnitude of the high density effects in the kinematical range of the future eA colliders. 

During the last years an intense activity in the area resulted in more sophisticated dipole proton cross sections, 
which had more theoretical constraints and which were able to give a better description of the more recent HERA 
and/or RHIC data p"8l - [24| . In what follows we will use the b-CGC model proposed in Ref. [l^, which improves 
the IIM model [1^ with the inclusion of the impact parameter dependence in the dipole proton cross sections. The 
parameters of this model were recently fitted to describe the current HERA data [lol • Following [l^ we have that 
the dipole-proton cross section is given by: 



,2\ _ 



= / d'b 



da dp 
IPb 



where 



dfj. 



dp 



d^b 



2J^P{x,r,b) = 2 X 



rQs < 2 

1 - exp [- a In^ (6rQs)] rQs>2 



(3) 



(4) 
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with Y = ln(l/a;) and k = x"{'ls)/x'{'ls), where % is the LO BFKL characteristic function. The coefScients a and 
b are determined uniquely from the condition that M'^{x^r) and its derivative with respect to rQg are continuous at 
rQs = 2. They are given by: 

and 6=^(1 -AAo)-^^ . (5) 



(l-A/-o)2 1n(l-AAo) 2 
In this model, the proton saturation scale Qs now depends on the impact parameter: 



Qs = Qs{x,b) 



r 



exp 



2ScGC / 



(6) 



The parameter i?cGC was adjusted to give a good description of the i-dependence of exclusive J/ip photoproduction. 
Moreover the factors Ao and 7^ were taken to be free. In this way a very good description of F2 data was obtained. 
The parameter set which is going to be used here is the one presented in the second line of Table II of [205: 7^ = 0.46, 
BcGC = 7.5 GeV-2, Wq = 0.558, xq = 1.84 x 10"^ and A = 0.119. 

In order to estimate the magnitude of the saturation effects in heavy quark production it is important to compare 
the CGC predictions with those associated to linear QCD dynamics. As a model for the linear regime we consider 
the leading logarithmic approximation for the dipole-target cross section [2^ , where adA is directly related to the 
nuclear gluon distribution xgA as follows 

adA{x, P) = -jr^CisXgAix, lO/r^) . (7) 

The use of this cross section in the formulas given below will produce results which we denote CT, from color 
transparency. In this limit we are disregarding multiple scatterings of the dipole with the nuclei and are assuming 
that the dipole interacts incoherently with the target. In what follows we consider two different models for the nuclear 
gluon distribution. In the first one we disregard the nuclear effects and assume that xgA{x,Q^) = A.xgpf{x,Q^), 
with xgj^ being the gluon distribution in the proton and given by the GRV98 parameterization [28]. We will refer to 
this model as CT. In the second model we take into account the nuclear effects in the nuclear gluon distribution as 
described by the EKS98 parameterization [2§|. In this case xgA{x,Q^) = A.Rg{x,Q'^).xgN{x,Q'^) with Rg given in 
[29| . We will call this model CT + Shad. In our calculations the charm quark mass is rric = 1.5 GeV and the bottom 
quark mass is mt, ~ 4.5 GeV. 



II. HEAVY QUARK PRODUCTION 



The electron-proton (ep) collider at HERA has opened up a new kinematic regime in the study of the deep structure 
of the proton and, in general, of hadronic interactions, which is characterized by small values of the Bjorken variable 
X — Q'^/s, where is the momentum transfer and ^/s is the center-of-mass energy. In this regime we expect that the 
usual collinea r ap proach jsoj be replaced by a more general factorization scheme, as for example the fc_L-factorization 
approach (3ll - l33j | or the quasi-multi-Regge-kinematics (QMRK) framework [s^l (For related studies see Refs. (35l - [37| ). 
Let us present a brief review of these distinct approaches. 

In the collinear factorization approach [13] all partons involved are assumed to be on mass shell, carrying only 
longitudinal momenta, and their transverse momenta are neglected in the QCD matrix elements. Moreover, the cross 
sections for the QCD subprocess are usually calculated in the leading order (LO), as well as in the next-to-leading 
order (NLO). In particular, the cross sections involving incoming hadrons are given, at all orders, by the convolution 
of intrinsically non-perturbative (but universal) quantities - the parton densities - with perturbatively calculable 
hard matrix elements, which are process dependent. The conventional gluon distribution g(a;, /i^), which drives the 
behavior of the observables at high energies, corresponds to the density of gluons in the proton having a longitudinal 
momentum fraction x at the factorization scale /i. This distribution satisfies the DGLAP evolution in /i^ and does not 
contain information about the transverse momenta k± of the gluon. On the other hand, in the large energy (small-x) 
limit, we have that the characteristic scale /i of the hard subprocess of parton scattering is much less than ^/s, but 
greater than the Aqcd parameter. In this limit, the effects of the finite transverse momenta of the incoming partons 
become important, and the factorization must be generalized, implying that the cross sections are now /c^-factorized 
into an off-shell partonic cross section and a /c^-unintegrated parton density function J^(x,k±), characterizing the 
fc^-factorization approach [sil - lssj . The function J- is obtained as a solution of the evolution equation associated to 
the dynamics that governs the QCD at high energies. A sizeable piece of the NLO and some of the NNLO corrections 
to the LO contributions on the collinear approach, related to the contribution of non-zero transverse momenta of the 
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incident partons, are already included in the LO contribution within the fc_L -factorization approach. Moreover, the 
coefficient functions and the splitting functions giving the collinear parton distributions are supplemented by all-order 
Qfs 111(1/2:) resummation at high energies [11]. A detailed comparison between the predictions of the collinear and 
/c^-factorization approaches for the heavy-quark photoproduction was performed in Refs. 0,11^, which we indicate 
for more details of these two approaches. 

In the last years, an alternative approach to calculated the heavy quark production at high energies was proposed 
considering the QMRK framework. It is based on an effective theory implemented with the non-Abelian gauge- 
invariant action obtained in Ref. (33 |. In this approach the initial-state i-partons are considered as Reggeons. In 
contrast to the fcj^-factorization approach, the QMRK approach uses gauge-invariant amplitudes and is based on a 
factorization hypothesis that is proven in the lea ding log arithmic approximation. The phenomenological implications 
of this approach were discussed in detail in Refs. |40l - t4^ . which demonstrated that the QMRK approach is a powerful 
tool for the theoretical description of the high energy processes. In particular, in the F| and Z?-meson spectra 
are successfully described using the QMRK approach. 

The heavy quark production can also be calculated using the color dipole approach (26l |. This formalism can be 
obtained from the /cj^-factorization approach after the Fourier transformation from the space of quark transverse 
momenta into the space of transverse coordinates (See e.g. (isj). It is important to emphasize that this equivalence is 
only valid in the leading logarithmic approximation, being violated if the exact gluon kinematics is considered A 
detailed discussion of the equivalence or not between the dipole and the QMRK approaches still is an open question 
(See, however, Refs. (45j). The main advantage to use the color dipole formalism, is that it gives a simple unified 
picture of inclusive and diffractive processes and the saturation effects can be easily implemented in this approach. 
It is important to emphasize that phenomenological models based on the Color Glass Condensate (See, e.g., [201) 
the solution of the running coupling BK equation (46, - ,48,] describe quite well the current experimental HERA data 
for inclusive and exclusive observables. 



A. Charm structure function 



In terms of virtual photon-target cross sections (Tt,l for the transversely and longitudinally polarized photons, the 
nuclear F2 structure function is given by 

F^ix,Q^) = -r-^<ytot{x,Q') 

47r Ctprn 



with (26|: 



crtot = (TT + (TL and arx = J (fr dz\'^T,Lir, z,Q'^)f adA{x,r), (8) 

where ^t,l is the light-cone wave function of the virtual photon and adA is the dipole nucleus cross section describing 
the interaction of the qq dipole with the nucleus target. In Eq. ([8]) r is the transverse separation of the qq pair and 
z is the photon momentum fraction carried by the quark (for details see e.g. Ref. 51]). The charm component of 
the nuclear structure function F2 (x, Q^) is obtained directly from Eq. ([8]) isolating the charm flavor. In Fig. [T] we 
show F2{x,Q'^) / A. As expected, in this kinematical domain it grows with and falls with increasing x. What is 
really remarkable is the difference between the bCGC and the linear CT models, which can reach a factor up to 4!. In 
previous estimates of this observable [1^, non- linear effects were found to be weaker. However, in [i^ the input was 
different (unintegrated gluon distribution instead of a dipole cross section) and the procedure adopted to estimate 
the purely linear contribution was to switch off the non-linear effects in the unintegrated gluon distribution of the 
proton. In some of our previous works (for example in Q) we adopted an analogous procedure and tried to make this 
separation switching off the non-linear component of the dipole-proton cross section. Although this method could 
give us some rough idea of the role played by some non-linear effects, we were missing part of them associated with 
the fusion of gluons belonging to different nucleons. Therefore we believe that the mentioned previous estimates have 
underestimated the importance of non-linear effects. 



B. Heavy quark spectrum 



As discussed before, heavy quark production has been very well studied over the last years both theoretically and 
experimentally. The elementary cross sections have been calculated in perturbative QCD up to next lo leading order 
and the parton densities have been extracted with the same degree of precision. Theory has been used to study heavy 
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quark production in ep collisions at HERA, in pp collisions at Tevatron and at RHIC, in pA and dA collisions at 
RHIC and in AA collisions at CERN-SPS and at RHIC. A recent and comprehensive survey of these advances can be 
found in [sO]- However, to the best of our knowledge, heavy quark production in eA has received almost no attention. 
This is probably related to the fact that eA data are old and until recently, there was no prospect of having high 
energy eA data. With the possible construction of a high energy EIC, updated estimates of heavy quark production 
are needed. We wish to address the subject from the perspective of saturation physics and thus the best option to 
obtain the production cross section, isolating shadowing and non-linear effects, is to use the dipole model. The dipole 
approach is very natural for the study of exclusive hidden charm and beauty electro and photo-production especially 
in the vector meson channel. As it was shown in j5l| the dipole formalism can be easily extended to open charm and 
beauty electro-production obtaining a quite successful description of the HERA data for the F| and D-meson spectra. 
Here, in order to calculate the differential heavy quark production cross section, da^'^ /d^pg, we have extended the 
approach of Ref. [5l|, which was originally developed for ep scattering, to electron-ion collisions with the Glauber- 
Gribov formalism. In this extension we implicitly assume that the factorization of the cross section verified for ep 
collisions remains valid in the nuclear case and make use of the dipole-nucleus cross section, which, in turn, depends 
on the dipole-nucleon cross section. For this last quantity we take the recent parametrization given by Eqs. ([3]) and 
The resulting cross section reads: 



GegUe 



da 



4gV(l-a)^ 



h 

4e 



Pcf + 



h 
2 
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(9) 



with 



h = J drrJo{pQr)Ko{er)adA{r) 
h = I drr^ JoipQr)Ki{er)adAir) 



= J drr Ji{pQr) Ki{er) adAir) 



(10) 



where Jo,i and Ko,i are Bessel functions, e ~ a{l — a)Q^ 



and adA is given by Eq. ([T]) or (O . 



To calculate the _D-meson production cross section we must let the charm quark fragment. Following |51| 
convolute the charm quark production cross section ^ with the nonperturbative fragmentation function: 



dcr{YA -> DX) f dp^ da da{j*A ^ cX) 



dz d'^pj=) J a d^p^ da 

where Dq{z*) is the well known Peterson fragmentation function given by 



Pd 



-Pc 



n{h) 



11 -i^ 



(11) 



(12) 



The fragmentation function gives the probability that the original charm quark with a momentum P fragments into a 
D-meson with momentum fraction z*P. There are more recent fragmentation functions but here we are only interested 
in checking if the differences between linear and non-linear dynamics are affected by fragmentation. For this purpose 
the Peterson fragmentation function is adequate. 

In Fig. [2] we show the transverse momentum spectrum of charm quarks. The main purpose of this figure is to show 
that the predictions of the linear physics (CT -I- Shad) differ from the total (i.e. bCGC) by a factor which increases 
with the energy W and goes from 1.5 {W — 100 GeV) to 4 [W = 1400 GeV). Moreover, this difference persists for a 
wide momentum window. At very large pr we enter the deep linear regime and expect that the two curves coincide. 

In Fig. |3]we show the transverse momentum spectrum of bottom quarks. As expected, we observe the same features 
of the charm distribution, except that now the non-linear effects are weaker. Nevertheless they are still noticeable. 
In Fig. |3]we show the dependence of the pr distribution at a fixed value px — ^ GeV^ for different energies. The 
upper and lower panels show the charm and bottom distributions respectively. Here again, we observe a remarkable 
strenght and persistence up to large virtualities of the differences between CT -I- Shad and bCGC. In Fig. [5] we show 
the transverse momentum spectrum of D mesons for three energies W = 200, 500, 1400 GeV and for two virtualities 
= 2 GeV^ (upper panels) and — 10 GeV^ (lower panels). As it can be seen, the differences between the curves 
CT, CT -|-Shad and bCGC are the same as before. 
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FIG. 2: (Color online) Transverse momentum charm spectrum for =2 GeV'^ and difTerent center-of-mass energies. 




III. CONCLUSIONS 

In this work we have updated the calculations presented in [sH and extended them to electron-ion collisions. We 
compared the predictions of a saturation model (bCGC) with the predictions made with a linear model. The main 
conclusion was that it seems quite possible to observe the non-linear effects both in the charm structure function 
F2'^{x, Q^) and in the px distributions of the heavy quarks. For the energies considered this difference is of a factor 
going from 1.5 to 4. As expected, the final state conversion of the heavy quarks into heavy mesons, performed 
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FIG. 4: (Color online) Dependence on the photon virtuality at — 4 GeV^ 



through the convolution of our px distributions with the Peterson fragmentation function, does not change the 
difference between "full" {— linear + non-linear) and linear predictions. In a future analysis of fragmentation we shall 
include the production of heavy mesons from light quarks. Although the fragmentation channel q ^ H (where q is 
a light quark and is a D or B meson) is disfavoured in comparison with Q H , the production of light quarks 
from the incoming photon is strongly enhanced by the photon wave function. One effect might compensate the other 
and, in the end, light quarks might play a significant role in heavy meson production [53 |. If this would be the case, 
non-linear effects would be even stronger. 

Our results suggest that heavy quark production in high energy eA collisions is a promising signature of saturation. 
Previous estimates of this same observable were not so positive, probably because they addressed ep or pp collisions 
as in Ref. (s^l or because the method employed to separate linear and non-linear effects was not very accurate. From 
our figures we can also conclude that non-linear dynamics, here as in several other contexts, leads to a depletion in 
the pt spectra, in contrast to some early estimates [s^l • 

A final comment is in order. As discussed in Section II, the heavy quark production at high energies can be 
calculated considering different approaches which are not equivalent in the full kinematic region. Consequently, a 
more detailed study of the saturarion effects using these different approaches is important in order to estimate the 
theoretical uncertainty of our predictions. It is postponed for a future publication. 
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FIG. 5: (Color online) Transverse momentum spectra of D mesons. 
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